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on Block Fading Channels 
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^ , Abstract 



The outage probability limit is a fundamental and achievable lower bound on the word error rate of 
coded communication systems affected by fading. This limit is mainly determined by two parameters: 
the diversity order and the coding gain. With linear precoding, full diversity on a block fading channel 
can be achieved without error-correcting code. However, the effect of precoding on the coding gain is 



' not well known, mainly due to the complicated expression of the outage probability. Using a geometric 

approach, this paper establishes simple upper bounds on the outage probability, the minimization of 

CO , which yields to precoding matrices that achieve very good performance. For discrete alphabets, it is 
> 

00 ' shown that the combination of constellation expansion and precoding is sufficient to closely approach the 

' minimum possible outage achieved by an i.i.d. Gaussian input distribution, thus essentially maximizing 



the coding gain. 



I. Introduction 

In many applications, using techniques such as frequency-hopping (GSM, EDGE), time- 
y ; interleaving (DVB-T), OFDMA (WiMax and LTE), H-ARQ with cross-packet coding [[T6]l. llSll. 

and cooperative communications [|29l , [l30l . [|20l . [|71, the channel can be modelled as a flat 
block fading (BE) channel where the fading gain is piecewise constant over the duration of 
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a transmitted packet. Due to motion of the transmitter, receiver or objects between the transmitter 
and receiver, the fading gains vary from one packet to the next and are considered unknown 
at the transmitter side. The fraction of codewords where decoding fails to wipe out all errors 
is referred to as the average word error rate (WER). When displayed on a log-scale versus the 
average signal-to-noise ratio (SNR) in decibel, the high-SNR slope of the WER is called the 
diversity ordeiQ. Since the diversity order determines how fast the error rate decreases with the 
SNR, it is then a key parameter of the communication system. 

By appropriately linearly precoding multidimensional signal constellations before transmitting 
each of its components on a different block of the block fading channel (denoted as component 
interleaving), the Singleton bound does not limit the maximal coding rate Rc, achieving full- 
diversity, any more [|T3l . [[T4l . [[T2l . In particular, full diversity can be achieved without error- 
correcting code. Therefore, it has been almost exclusively studied for uncoded schemes (see [HI, 
[El and references therein, e.g. [O), except for few papers investigating outer coded transmission 
schemes for MIMO (e.g. [[T3l . |fT9ll ) and recent work on pairing subchannels when channel 
state information at the transmitter (CSIT) is available [[23l . [[24|. However, besides yielding 
full diversity with a higher rate, linear precoding can also improve the coding gain of coded 
systems. So the effect of linear precoding on the diversity order is well understood, but there is 
no sufficient insight in the effect of linear precoding on coding gain of coded systems. 

Before studying the optimization of the WER for practical schemes with linear precoding, it 
is important to understand the influence of linear precoding on the performance limits of the 
communication channel. The outage probability limit is an achievable lower bound on the WER 
of coded systems in the limit of large block length [0, [[25l and is given by ll32l Sec. 5.4], 

Pout(7,i",^) = Pr(/(«,7,^)<i?), 

where /(q;,7,P) is the instantaneous mutual information between transmitted symbols and 
received symbols as a function of the set of fading gains a. observed during the transmission of 
one codeword, the average SNR 7, and the precoding matrix P. By choosing a well designed 
precoding matrix P, the outage probability can be minimized. But only a brute force optimization 

'We assume fading gain distributions whiere the WER can be expressed as gSNR^'' where g and 5 are constants; and 5 is 
the diversity order. 
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can minimize the outage probability as a closed form expression of the outage probability is 
not available. Such an optimization is often intractable when the number of fading gains per 
codeword is larger than two and/or large constellations (e.g. 16 points) are used. A simple 
approach could be designing P such that the mean of / (ct, 7, P) over the fading distribution 
is maximized, in the hope that the area under the left tail of its probability density function 
(pdf) would be minimized. However, the ergodic mutual information Eq[/ (o:, 7, P)] contains 
no information on the diversity order, and, due to the limited spectral efficiency of finite discrete 
input alphabets, for increasing 7 rapidly converges to a maximum that does not depend on the 
precoding matrix. Therefore, the maximization of the ergodic mutual information fails to provide 
the optimum precoding matrix [|9l. 

In this paper, we first study the effect of linear precoding of discrete input alphabets on the 
outage probability^. A unitary B x B precoding matrix has B'^ parameters that are free to choose 
(denoted as degrees of freedom), so that the optimization of the precoding matrix is multivariate. 
In a brute force optimization, Monte Carlo simulations are required to take into account the 
distribution of all fading gains when computing the outage probability. Here, the analysis uses 
a geometric approach that leads to upper bounds on the outage probability, which are easier to 
optimize, because it is no longer necessary to perform a Monte Carlo simulation based on the 
fading gain distribution. Therefore, with a minimal computational effort, the precoding matrix 
minimizing the upper bound on the outage probability limit can be determined. These results 
serve as a basis for the design of practical coded systems with linear precoding [|6l, which is 
discussed in the final section before conclusions. 

II. System model 

At the transmitter output, a packet is represented as a real- valued or complex column vector 
X = • • • ) of dimension N, consisting of B blocks that each contain N/B sym- 

bols, where (.)^ indicates transposition; the b-th block of the packet is x(6) = [xib)i, . . . , 
with E[|x(6)„p] = 1. The channel is memoryless with additive white Gaussian noise and 
multiplicative real-valued fading. The fading coefficients are only known at the decoder side. 

^Note that for i.i.d. Gaussian input distributions, it was found that a power allocation (coixesponding to a multiplication of 
the symbol vector with a non-unitary matrix) over the different blocks minimized the outage probability |311 Sec. 5.1], 1171 . 
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The received signal vector corresponding to the transmitted block x{b) is 



fJ-ib) = abX{b) + i^{b), 6=1,... 



B. 



(1) 



The fading coefficients {ab, b = 1, . . . , B} are i.i.d. In the numerical results, we consider the 
fading coefficients to be Rayleigh distributed, with ]E[a^] = 1, but the analysis in this paper does 
not depend on the fading distribution. The noise vector uj(b) consists of N/B independent noise 
samples which are complex Gaussian distributed, uj{b)n ~ CA/'(0, 2cr^). The average signal-to- 
noise ratio is 'J = 7r2- 

The transmitted vector x is obtained from the information bits through a sequence of opera- 
tions. Assuming a binary encoder with coding rate Rc, a packet of K information bits is encoded 
into K/Rc coded bits. The binary codeword is split into K/{mRc) strings each containing m 
bits. In a standard coded communications system, the components of the transmitted vector 
x{b) are obtained by directly mapping each string of m coded bits to one of M = 2™ points 
belonging to a 1 -dimensional real or complex space; the corresponding spectral efficiency R 
in bits per channel use (bpcu) is given hy R = mRc. When using precoding combined with 
component interleaving, each string of m bits is mapped to one of M = 2™ points belonging to 
a B-dimensional real or complex space; the corresponding B-dimensional M-point constellation 
fi^ is denoted M-TZ^ or M-C^, respectively. Denoting as z(n) = [2(72)1, ^(n)^]^ the B- 
dimensional vector that results from mapping the n-th string of m coded bits, the linear precoding 
involves the computation 



where P is a non-singular precoding matrix of dimension B x B. The precoder output vectors 
x(n) = . . . ,x{n)B] belong to a B-dimensional M-point constellation fi^ which results 

from a linear transformation (through P) of Q^. Finally, component interleaving is applied: 
the n-th component of the b-th block of the transmitted vector x equals the b-th component 
of the n-th precoder output vector, i.e., xib)n = x{n)b (Fig. [l]). Hence, the B components of 
x(n) experience independent fading when transmitted over the BF channel. Taking into account 
that K information bits are transformed into a transmitted vector x with = KB/{mRc) 
components, the overall spectral efficiency is i? = R^.^ bpcu. Note that there are several ways 
to achieve a given spectral efficiency R. For example, R = 0.9 bpcu for i? = 2 can be achieved 



x(n) = Pz(n), n = 1, . . . , 



N 

'b 



(2) 
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Fig. 1. The coded bits are mapped to multidimensional symbols z(n) in a constellation flz- Each symbol z(n) is transformed 
to x(7i), where the different components are placed in the corresponding blocks in x- 



by a coded communication system with precoding and component interleaving (R = Rc^) by 
choosing m = 3 and Rc = 0.6, whereas a standard communication system (R = Rem) achieves 
R = 0.9 bpcu for m = 1 and R^. = 0.9. Also other combinations of m and Rc are possible. 
Note that a standard communication system with R = rn'R^ can be viewed as a special case of 
a multidimensional modulation, where is a Cartesian product of B identical constellations 
belonging to 2™' -7^^ or 2™'-C^ and P = I, where / is the identity matrix, so that il^ = ^z', 
these B -dimensional constellations contain M = 2"^'^ points. 
We can reformulate Eq. ([T) in terms of x(n) as 

y(n) = a ■ x(ra) + w(n), n = l,...,—, (3) 

where y{n)h = /i(fe)n, {oc)b = ctb, w{n)b = co{b)n and ex ■ x(?t,) denotes component-wise 
multiplication: (ck ■ x(n))b = abx{n)b. 

As precoding allows to increase the spectral efficiency associated with full diversity, this 
technique has been extensively studied in previous works, but almost exclusively for uncoded 
schemes [|2l, Q. Here, we study this system for coded schemes and we choose the precoding 
matrix P and the constellation Vt^ minimizing the outage probability. 

A precoding matrix P that is unitary is a natural choice because it does not decrease the 
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capacity of a Gaussian channel. In this paper, we restrict our study to real-valued precoding 
matrices, hence P is orthogonal. When i? = 2, P is a rotation matrix where the rotation angle 
6 is the degree of freedom: 

cos(6') — sin(6') 
sin(6') cos(6') 



(4) 



However, rotation matrices are difficult to construct for higher dimensions. In Sec. UlIl it will 
be shown that for B > 2 it is sufficient to consider orthogonal circulant precoding matrices. We 
denote its first row as (po, • • • The second row is a cyclic shift to the right of the first 

row, and so on. Because the columns of the B x B Fourier matrix F are the eigenvectors of 
any circulant matrix, we can construct P as follows: 



FAF 



H 



(5) 



where (F)m,n = -^exp (■ 



-2j7rmn ^ 



^ m, G {0, . . . , i? — 1}, and A is a diagonal matrix containing 
the eigenvalues of P. The condition for P being orthogonal is A^A = Ib, or the B eigenvalues 



of P must have a squared magnitude of 1. It is easy to find that 

B-l 



A. 



Ef —j2Txnl\ 
1=0 ^ ^ 



Now, it follows that 



Pi 



E ( j27Tml\ 
^ A„exp(^-^J 

m=0 ^ ^ 



(6) 



(7) 



As the eigenvalues must have a magnitude of 1, we have A„ = exp(j6'„). In order to obtain 
a real- valued P, we take Aq real- valued (i.e., Aq = 1 or Aq = —1) and Xb-u = (An)* (i-e., 
Ob^u = —dn) for n = 1, . . . , P — 1. When P is even, this implies \b/2 = 1 or \b/2 = —1- For 
P > 2, P is determined by [(P — 1)/2)J continuous parameters that can be optimized. 

Note that for P = 3, P constructed as above, with Aq = 1 and Ai = exp(j6'i), corresponds 
to a 3-dimensional rotation with angle 9i around the fixed axis ^(1,1,1), 



P 



1 
3 



l-k + ^/Sl l + 2k l-k-y/3l 
l-k-Vsi l-k + V3l l + 2k 



(8) 
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(a) (b) 



Fig. 2. Displaying the rotation at the transmitter |(a)| and at the receiver (without noise) [(b)] for B = 2. The empty (filled) 
circles represent (ilx)- The components of t(n) are obtained by scaling the components x(n) by their respective fading 
gain (here Q2 < Qi), as expressed by the component-wise multiplication a ■ x(n). The crosses on the coordinate axes are the 
transmitted |(a)| and received |(b)| vector components, respectively. 



where k = cos(6'i) and / = sin(6'i). 

Fig. [2] illustrates the effect of a rotation for i? = 2 when a 4-7^^ constellation is used as fi^. 
The transmitted components are affected by their corresponding fading gain, which is expressed 
by the component- wise multiplication t(n) = a ■ x(n), which is shown at the right side in Fig. 
m We say that t{n) belongs to the faded constellation The point t(n) = [t{n)i, . . . ,t{n)B] 
is shown for a particular fading point a. When ai ^ a2, the constellation fi^ can be interpreted 
as a distorted QPSK constellation (i.e., a constellation in which both components do not have 
the same magnitude). 

Considering t(n), an equivalent channel model can be formulated: 

y(n) = t(n) + w(n), n = l,...,— , (9) 

which yields more insight and will be useful in the proofs of propositions of this paper. 
This system model represents a Gaussian vector channel with input t(n). This means that for 
a particular fading point, the block fading channel can be interpreted as a virtual^ Gaussian 
channel, with a discrete input alphabet fi^. 

^We use the term virtual because the fading gains of the actual channel are incorporated in the constellation Qt- 
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III. Analysis of Outage Probability in the fading space 

For the remainder of the paper, we will drop the index n in the vectors z(n), x(n), t(n), 
y(n) and w{n), as the time index is not important when considering mutual information. We 
write random variables using upper case letters corresponding to the lower case letters used 
for their realizations. The mutual information / (a, 7, P) at a certain fading point ct between 
the transmitted B-dimensional symbol x (uniformly distributed over il^) and the corresponding 
received vector y is given by 

/(a,7,P) = i/(X;Y|a,7), (10) 

where the scaling factor is added because the B blocks in the channel timeshare a time-interval 
E Section 9.4], [32, Section 5.4.4]. The mutual information /(X; Y|q:,7) is [HI, ifTOll 

/(X;Y|a,7)=m-2-^'"^Ey|, 

where (i^(v,u) = J2b=i 1^^ ~ "^bf- However, more insight can be gained when considering t. 
From Eq. it is clear that for a certain fading point, the mutual information of this virtual 
channel, with input ilt, is the same as the mutual information of the actual channel, with input 

/(a,7,P) = -^/(T;Y|a,7). (12) 

Hence, the fading point o: maximizing (minimizing) the mutual information /(X; Y|q;, 7) corre- 
sponds to the fading point that distorts constellation fi^ in the best (worst) way at the input of a 
Gaussian vector channel. This interpretation allows to exploit the many results from literature on 
the mutual information of Gaussian channels. Therefore, Eq. (fT2)) will be useful in the following. 

The outage probability is the probability that the instantaneous mutual information is less than 
the transmitted information bitrate R ll32l section 5.4], 

Pout(7,^,i?) =Pr(/(«,7,i') < R)- (13) 
Our main goal is to find the precoding matrix P that minimizes the outage probability Pout(7, P, R), 

Popt = argmin Pout (7, P, R)- 
p 



logs I 



2a2 



(11) 



January 18, 2013 



DRAFT 



9 




(a) Examples of an outage boundary when — 4- (b) The points {0^,0, b = 1,...,B} and are 
TZ^ and the rotation angle is ^ = 0, 10, 27 degrees. shown for B = 2. 
The region K is coloured red for 9 — 27 degrees. 



Fig. 3. The outage boundary limits the outage region K in the fading space which corresponds to an information theoretic 
outage event. In the Figure, the information rate is _R = 0.9 bpcu and the average SNR is fixed to 7 = 8 dB. 



A closed form expression for J(q:,7,P) does not exist and Eq. (fTTj) is difficult to analyze 
because of the presence of the mathematical expectation. Therefore, in the following two sections, 
we will develop simple bounds on the outage probability and will approximate Popt by the optimal 
precoding matrix that minimizes the bounds on the outage probability. 

But first we will introduce a framework in this section that allows to gain insight on the 
meaning of the outage probability. The considered framework is the fading space JH, which is 
the B-dimensional Euclidean space, , of the real-valued positive fading gains. The outage 
probability equals the probability that the fading point ct belongs to a region Vo, which is such 



that / (a, -f,e) < R for all cx in K (Fig. 3(a)): 



Pout(7,^,^) = / p{a)da, (14) 
where p{oc) is the joint pdf of the fading gains ai, . . . , a^. We say that the region K is limited 



by an outage boundary Po(7, P, R) (see Fig. |3(a)| for B = 2), defined by 



/(a,7,P) = P, Va G Po(7,^,^)- 
In general, we can consider a certain boundary in the fading space limiting a certain region V. 
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This boundary is denoted by B(y). The outage boundary i?o(7, P, R) corresponds to the region 
Vo, as defined in (fT4l) . so that the outage boundary is -Bo(7, -P, R) = -B(K)- 

Definition 1: We define a^^o by the magnitude of the intersection between the outage boundary 
and the axis ai,. More precisely, / [a.\a-=oA^b,ai=ato^'y^ R) = R- ^y convention, o = +00 if 
the axis is an asymptote for the outage boundary. It is possible that o does not exist. 

Definition 2: We define as the value of the components of the intersection between the 
outage boundary and the line ai = . . . = as (also known as the ergodic line). More precisely, 

/ (Q:|a,=ae,iG[l,...,i?],7,^) = R- 

The defined points are illustrated in Fig. |3(b)| for 5 = 2. In the remainder of the paper, we 
denote the points a\a,=o,i=^b,ab=ai,, by a^^o and a\a,=a„ie[^,...,B] by cte- 

Proposition 1: On a BF channel with B = 2, with a discrete input alphabet and with linear 
precoding, the magnitudes {a;;,o, b = 1,2} are equal if the constellation is invariant under a 
rotation of 90 degrees and the precoding matrix is orthogonal. 

On a BF channel with B > 2, with a discrete input alphabet and with linear precoding, the 
magnitudes {ab,o, b = 1, . . . , B} are equal if the constellation is invariant under a cyclic shift 
of the components of the points of the constellation and the precoding matrix is an orthogonal 
circulant matrix. 

Proof: See Appendix lAl ■ 
Notice that the condition of Prop. 1 is a sufficient condition and not a necessary condition. In the 
remainder of this paper, it is assumed that the constellation used at the transmitter fulfils Prop. 
[TJ The magnitudes {ab,o, b = 1, . . . , B} will then simply be denoted by ao- This also means that 
the projection of fi^ on either coordinate axes yields the same set of points, which we denote 
by Sp, where p stands for projection. 

Note that a multidimensional constellation that fulfils Prop. \T\ (i.e., its projection on each 
coordinate axis yields the same set of points, see Appendix lA]) has an interesting property. 
For these constellations, the function / {Xf,; Yf,\ab = a, 7), which is the mutual information of a 
point-to-point channel with fading coefficient ab = a, average SNR 7 and with discrete input 
Xb, does not depend on b. As a consequence, we will represent this mutual information by 

Definition 3: Boundary Bi = B{Vi) is said to outer bound outage boundary B2 = B(y2) if 
V2 C Vi. The opposite holds for inner bounding. 
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The outage boundary in the fading space has a simple but interesting property: if an outage 
boundary outer bounds (inner bounds) another outage boundary, than its corresponding out- 
age probability is larger (smaller) (see Eq. dHJ). For example, consider the input alphabet 
Z ~ A/'(0,/), which we denote as an i.i.d. Gaussian input alphabeu where is a column 
vector of zeros and / is the identity matrix. We denote the outage boundary (outage region) 
corresponding to i.i.d. Gaussian inputs and discrete input alphabets, by 5o(Gauss) (K(Gauss)) 
and i?o(discrete) (K(discrete)), respectively. The boundary i?o(Gauss) inner bounds i?o(discrete) 
jSl. Therefore, the outage probability corresponding to i.i.d. Gaussian inputs is a lower bound on 
the outage probability corresponding to a discrete input alphabet. Consequently, by minimizing 
the outage probability corresponding to a discrete input alphabet, i?o(discrete) can at most 
approach i?o(Gauss) (see Fig. H]). 



Gaussian distributed continuous input 
Discrete uniform distributed input 




Fig. 4. The outage boundary _Bo(discrete) is inner bounded by _Bo(Gauss). Optimizing the precoding matrix can at most make 
_Bo(discrete) approach Bo(Gauss) as illustrated by the arrows. 



In the following two sections, we will determine boundaries with simple shapes outer (inner) 



''We mainly refer to real-valued i.i.d. Gaussian input alphabets in this paper. In Sec. IVII-A2I we extend this to complex i.i.d. 
Gaussian input alphabets. 
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bounding i?o (discrete), which are then much easier to optimize. For the outer boundary, this can 
be done by determining a surface in the fading space, U (a) = 0, satisfying 

I (cK, 7, P) > R, for all a satisfying U{ct) = 0. (15) 

For the inner boundary, the greater than or equal sign is replaced by a smaller than or equal 
sign. For example, for B = 2, we will prove that a circular arc touching i?o(discrete) at ai = ao 
(on the horizontal axis) and a2 = Oo (on the vertical axis) satisfies Eq. (fTSi) . 



IV. Bounds on the outage probability without linear precoding 

As a first step, we will establish upper and lower bounds on the outage probability of a 
communication system without linear precoding. This will set the stage for the bounds with 
linear precoding in Sec. |V] 

We will prove that the outage region is outer bounded by a i?-hypersphere touching the outage 
boundary on the axes of the fading space, hence with radius ao- A 5-hypersphere U (a.) = is 
a generalization of a sphere to B dimensions, 

B 



2 2 



6=1 

Note that this is only possible for constellations fulfilling the conditions of Prop. \T\ For other 
constellations, the outer bounding _B-hypersphere will have a radius of maxftg[i ...^b] Ofe.o and will 
therefore be less tight. 

Next, we will prove that the outage region is inner bounded by a _B-hypersphere touching the 
outage boundary at a.^. 

Lemma 1: On a BF channel with a discrete multidimensional input alphabet X, the mutual 
information / (a, 7, P) is upper bounded as follows 

1 ^ 

na,^,P)<-J2^sA^h,P)- (16) 

6=1 

Proof: See Appendix |B] ■ 
A special case is the case without precoding (P = I) so that ^Ij. = ^1^. If in that case is the 
Cartesian product of il^^, . . . , ^Izg u> = i<Sp)^, then all variables Xt, = are independent 

^We denote the projection of on the 6-th coordinate axis as Q^^ - For constellations fulfilling Prop. [T] fi^j, = Sp. 
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for all b, so that equality holds in Eq. (fT6l) . E.g., fi^ = 4-7^^, as shown in Fig. |2l is the Cartesian 
product of two BPSK constellations, denoted as BPSKxBPSK. 

For an i.i.d. Gaussian input alphabet, the instantaneous mutual information with linear pre- 
coding is the same as the instantaneous mutual information without linear precoding, because 
Pz = X ~ A/'(0, /) when P is orthogonal, so that precoding does not change the distribution 
of the input vector. Therefore, the instantaneous mutual information for an i.i.d. Gaussian input 
alphabet is 



Proposition 2: The outage region ^^(Gauss) of a BE channel with an i.i.d. Gaussian input 
alphabet is outer bounded by the B -hyper sphere + a| + . . . + a| = and inner bounded 

by the i?-hypersphere + a| + . . . + a| = 

Proof: In Appendix O it is proved that for a on a hypersurface: 



• ^(o:,7) Ui=...=aB > ^(",7) 

because the mutual information is a sum of concave functions of the instantaneous SNRs 70^ 
(Eq. (flTl) ). Calculating cto and a^, according to Definitions [H & [2l yields the radii of both 



The inner boundary touches the outage boundary in the point CKe, so that it does not depend on 
B (the entire codeword is affected by the same fading gain a^). The outer boundary touches the 
outage boundary in the points cxbp, V b; its dependence on B follows from the fact that B — 1 
is equal to the number of erased channels. 

Proposition 3: On a BE channel with a discrete input alphabet that is a Cartesian product 
of one-dimensional constellations, the outage region ^(discrete) is outer bounded by the B- 
hypersphere al + + . . . + a% = al touching it at the axes of the fading space at cxb o, V b. 
Also, K(discrete) is inner bounded by the 5-hypersphere af + al + . . . + a% = Bal touching 
it at dis- 
proof: Lemma[i|proved that the mutual information is upper bounded by Xlfc=i ^Sp («b7) P)^ 
where this upper bound coincides with the exact expression in the case that Vt^ is a Cartesian 
product. Using the relation between the mutual information /(SNR) = Is^ {olIi^ P) ^^'^ ^he 
minimum mean-square error (MMSE) in estimating the input symbol, E Sp, given the output 
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(17) 



6=1 



5-hyperspheres. 



14 



symbol n ESI, 



d 



/(SNR) = -MMSE(SNR), 
2 



(18) 



d SNR 



it is easily proved that I^^ (ctl'y, P) is a concave function (the second derivative is negative) of 
the instantaneous SNR 70^, because the MMSE is a decreasing function of the SNR. Therefore, 



Notice that the techniques of the proofs of Props. |2] and |3] cannot be used for discrete input 
alphabets that are not a Cartesian product of one-dimensional constellations, because in that 



mutual information. However, this case is merely a particular case (P = /) of a preceded discrete 
input alphabet, which is covered in the next section. 



Propositions [2I and |3] mainly state that the outage region for a channel with an i.i.d. Gaussian 
input alphabet or a discrete input alphabet without precoding is outer (inner) bounded by a 
5-hypersphere touching it at a^oj V h (cKg). We conjecture that this property still holds for a 
communication system with a discrete alphabet with linear precoding at the input of the channel. 
First, we will give new detailed proofs for low and high instantaneous SNR of this property. 
Then, a more intuitive explanation will be given to provide more insight. 

Proposition 4: On a BF channel at low instantaneous SNR, with a discrete input alphabet 
and with linear precoding, the outage region K(discrete) is outer (inner) bounded by the B- 
hypersphere \ \ . . . \ a\ = a\ (a\ + a| + . . . + «| = fia^). 

Proof: See Appendix ■ 

Proposition 5: On a BF channel at high instantaneous SNR, with a discrete input alphabet 
and with linear precoding, the outage region \4(discrete) is outer (inner) bounded by the B- 
hypersphere af + + ■ ■ ■ + a% = al (a^ + 0^2 + . . . + a| = Ba^). 

Proof: See Appendix |El ■ 

The outer and inner bound touch the outage boundary at the axes of the fading space and on 
the ergodic line, respectively. The outer (inner) region corresponds to an upper (lower) bound on 
the outage probability. Minimizing this upper (lower) bound is simply achieved by minimizing 

fto (tte)- 



the proof is the same as for Prop. |2l 



case the upper bound 



^6=1 ^Sp i^bl^ ^o^^ '^ot coincide with the exact expression of the 



V. Bounds on the outage probability with linear precoding 
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To get more insight, we give an illustration of Props. |4]and |5] We consider the 4-7^^ constel- 
lation from Fig. |2(a)[ with 9 = n/6. Further, we take ai = aoCos(A) and 0^2 = aoSin(A), 
i.e., |a| = tto, so that a is on the outer boundary of the outage boundary. Fig. |5] shows the 
corresponding mutual information /(a, 7,^) as a function of 7 for various values of A. We 
observe that the mutual information increases when A increases from to 7r/4. The projections of 
the constellation fi^ on the horizontal and vertical coordinate axes have variances af = al cos^(A) 
and = a^sin^(A) respectively, whereas the total variance of ^It equals a^, irrespective of A. 
When A = 7r/4, is equivalent to a common QPSK constellation: the variances of both 

2 

projections equal and the corresponding mutual information is maximum. When A decreases 
from 7r/4 to 0, the difference between the two variances increases, Vlt becomes an increasingly 
more distorted QPSK constellation, and the mutual information decreases. At A = we get 
«! = ao and ^2 = 0; fit (which is now a scaling of iSp) reduces to a 4-PAM constellation 
and the corresponding mutual information is minimum. This illustrates the general principle that 
performance is optimized when the transmit power is equally split over the available identical 
channels and the performance is worst when the transmit power is completely used for only 
one channel, which is exactly what is claimed in Props. |4]and|5]for high and low instantaneous 
SNR. More specifically, for a on the hypersurface of a hypersphere, it is proved in Props. 5] 
and [5] that the mutual information is the smallest at and the largest at a.^. 

Prop. [6] states the condition under which the precoded constellations achieve full diversity for 
a rate R. This is proved in ffT2l . but based on Prop. |5l a new simple proof can be given. 

Proposition 6: For any coding rate < Rc < 1, the outage probability of a Rayleigh dis- 
tributed block fading channel with B blocks, with a discrete input alphabet and linear precoding 
exhibits full diversity, i.e., Pouti'J,P,R) oc 1/7^ at large 7, if 2^^ does not exceed the number 
of points contained in the projection of the precoded input constellation fi^ on a coordinate axis. 
Proof: See Appendix E ■ 

It should be noted that the proof of Prop. |6] assumes the existence of I^^{BR, P); hence, the 
number of points in Sp must not be less than 2^^. We consider the 4-7^^ constellation from 
Fig. |2(a)| to illustrate the constraint on the existence of I^^{BR, P). For 4-7?.^ from Fig. |2(a)| 
Sp contains only 2 points when the rotation angle ^ is a multiple of n/2, 3 points when the 
rotation angle 6^ is a multiple of 7r/4 but not of n/2, and 4 points otherwise, yielding R = 0.5, 
R = 3/4: and R = 1, respectively, as maximal rates corresponding to full diversity. 



January 18, 2013 



DRAFT 



16 




2 



4-QAM 
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Fig. 5. The mutual information is monotonically non-increasing as the constellation becomes more distorted. The upper and 
lower curve correspond to A = 7r/4 and A = respectively. 



VI. Optimization of the outage probability of precoded constellations 

In the previous section, we established that, for high and low instantaneous SNR, the outage re- 
gion of block fading channels with precoded constellations is outer bounded by a 5-hypersphere 
with center in the origin and touching the outage boundary on the axes of the fading space. 
This hypersphere corresponds to an upper bound on the outage probability of the channel (see 
the paragraph after Def. |3]). Instead of minimizing the actual outage probability, it is easier to 
minimize the upper bound on the outage probability. This optimization will allow the actual 
outage probability to closely approach a lower bound on the outage probability, i.e., the outage 
probability corresponding to an i.i.d. Gaussian input alphabet, as we will see in the numerical 
results. 

The 5-hypersphere is completely determined by one variable, its radius tto- We denote the 
i? -hypersphere outer bounding the outage region by Kp(tto) and the corresponding upper bound 
on the outage probability by Pup(tto): 




(19) 
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^2 



(a) The transmitted symbols are the components xi, (b) The received symbols t(, are the faded compo- 
of the constellation Q.^- nents abX^j. Collecting the different components ti,, 

the faded constellation Q,t is constructed, which is 

expressed hy Q.t = ol ■ Q.x- 

Fig. 6. The transmitted {^x) and faded constellation (ilt) are shown for the case of the transmission of real-valued symbols, 
B — 2 and ^1^ — For fit = a • Qx, the fading point (qi.o, 0), the intersection between the outage boundary and the axis 

02 ~ 0, is used. 



From Eq. (fT9l ). it is clear that the region Kp(«o) has to be made as small as possible to minimize 
Pup(tto)- Therefore, the optimization target is to minimize the radius ao- 

Because Isp{al'j,P) = BR, the minimization of ao (and, therefore, the minimization of the 
upper bound on the outage probability) is achieved by selecting the constellation Sp requiring 
the least energy to achieve a rate BR on a Gaussian channel. This involves a proper selection 
of both the constellation ^1^ and the precoding matrix P. Note that this optimization is much 
simpler than the direct minimization of the outage probability as given by Eqs. (fT3l) and (fTOl) . 
because the outage probability is hard to evaluate, especially when the number of fading gains 
and constellation points is large. Furthermore, no insight is gained by the latter approach, so 
that it would not be clear which constellation fi^ should be taken. 

Also, the outage region of block fading channels with precoded constellations is inner bounded 
by a B-hypersphere with center in the origin and touching the outage boundary in the point 
CKg. Hence, determines the radius of this fi-hypersphere, \fBa^. In this point, the faded 



constellation fif is balanced, as in Fig. |6(a)[ We will denote this balanced constellation Vtt = 
CKg ■ by Se- Hence, the hypersphere corresponds to a lower bound on the outage probability 
which is minimized by optimizing the mutual information of S^. 

In this section, we will show that the radius of the outer region, cto, and the radius of the 
inner region, ^/Ba^, can be minimized by combining a simple optimization of the precoding 
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matrix P with a constellation expansion. 

A. Optimization of the precoding matrix 

Precoding corresponds to performing a unitary transformation of the input vector x. As a 
unitary transformation preserves distance, it follows from Eq. (fTTI) that the value is insensitive 
to orthogonal transformations. Hence, the selection of the precoding matrix P affects only the 
radius of the outer region a^- Let us denote by O the set of parameters from P over which we 
will minimize Oo- For B = 2 and B = 3, the only degree of freedom is the rotation angle (see 
dl]) and dH])). For B > 3, more degrees of freedom can be exploited to minimize cto- For the 
numerical results, we restrict ourselves to < 3. 

The mutual information of iSp can be rewritten as Ig^ya^'^, O), which yields = ^ . 

Changing the value of O (e.g. the rotation angle 6 iox B = 2) will change the distances between 
the points in and so change its mutual information. For a fixed spectral efficiency R and 
fixed average SNR 7, minimizing the radius yields the optimization criterion 

Oopt = argmin O). 

o ^ 

The optimization is performed by means of a simulation, due to the lack of closed form expres- 
sions of the mutual information. Because the constellation is one-dimensional, the computational 
effort is minimal. We apply this optimization for different scenarios in Sec. IVIII 

B. Constellation expansion 

As the number of information bits per channel use is i? = mRc/B, there are different 
combinations of m and yielding the same R. Taking into account that R^ < 1, the minimum 
value of m equals \BR], with a corresponding coding rate Rc = jffj- 

The number of points in the constellation is \^lt\. Increasing the constellation size of will 
render a constellation fi^ with more points, both for and Sp. This higher order constellation 
may need less energy to achieve the same rate, both for the balanced case (optimization of 
fte) as the distorted case (optimization of ao). However, the decoding complexity increases as 
well as the complexity of optimization, so that there is a trade-off between performance and 
complexity. The higher the constellation size, the smaller the horizontal SNR-gap between the 
outage probabilities corresponding to a precoded discrete input alphabet and i.i.d. Gaussian input 
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alphabet. However, the improvement in performance becomes smaller and smaller, as illustrated 
in Sec. IWl 

VII. Numerical results 

A. Numerical results for B = 2 

When B = 2, O = 9, and the optimization criterion for the upper bound on the outage 
probability is to find 9 so that I^^{BR,9) is minimized. Next, a constellation expansion is 
performed to further minimize the upper bound as well as the lower bound on the outage 
probability. 

1) Real constellations: Assume that a transmission rate R = 0.9 bpcu is aimed. First, we 
consider the optimization of the rotation angle 9, see Fig. Ul On the left y-axis, we show the 
instantaneous SNR per symbol, 7s = a'^'j, so that Isfils,9) = BR. The minimum SNR per 
symbol 7s that is needed to transmit R = 0.9 bpcu for 7 = 8 dB is achieved by an i.i.d. 
Gaussian input alphabet: 

7s = 7; • 



This fundamental minimum can be approached when using a precoded discrete input fi^ = 4- 
TZ^ (Rc = 0.9) with rotation angle 9 = 27 degrees. Now, we apply a constellation expansion to 
further reduce ao (see Fig. |7]) and (see Fig. |8(a)[ ). For example, 7s for the rotated constellation 
= S-'R? {Rc = 0.6) approaches the theoretical minimum very closely for rotation angles 
within [0 . . . 9] degrees. An expansion to il^ = 16-7^^ (Rc = 0.45 and 6'opt G [35,45] degrees) 
only slightly improves the performance. The optimization of 7s decreases the volume of the 
region Vo, which is illustrated in Fig. |3(a)| for fi^, = 4-7?.^. Fig. |8(a)| illustrates that constellation 
expansion is sufficient to reduce the value a^, which is very close to the theoretical minimum. 
Therefore, the constellation is not further shaped to minimize the lower bound of the outage 
probability. 

The information theoretic approach used in this paper does not lead to the same optimized 
rotations as in the case of algebraic constructions of uncoded constellations. In [3 J and [1], 
multidimensional rotations have been optimized for uncoded infinite constellations transmitted 
on ergodic fading channels. As a simple illustration, we show in Fig. |7] the minimum product 
distance dp ^nm Il32]| of the uncoded S-TZ^ versus the rotation angle. The optimum and the profile 
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for Rotated 1 6-3!" + 

Ys for Rotated 8-3t^ x 

Ys for Rotated 4-3t^ * 

Ys for Gaussian 



Fig. 7. The optimization of thie radius of the outer region is shown for fl^ — R — 0.9 bpcu and 7 = 8 dB. The y-axis at 

the left denotes the instantaneous SNR per symbol, 7., = a'^'y, and the right y-axis denotes the minimum product distance. The 
thick black line without markers represents the fundamental minimum SNR per symbol, 7,, that is needed to transmit i? = 0.9 
bpcu for 7 = 8 dB, i.e., when using an i.i.d. Gaussian input alphabet. The effect of constellation expansion on the radius Oo 
is shown by optimization of f2z — and 0,^ = IQ-TZ^. The profile of the minimum product distance 

= 8-7^^ 



dp.min, is shown for 




Fig. 8. The outage boundaries (left) and outage probabilities (right) of Q,z = 4-7?,^, Qz ~ and f2z = 16-7?,^ with and 

without optimized rotation angle are shown. The spectral efficiency is i? = 0.9 bpcu and 7 = 8 dB (left). 



of (ip,min and those of 7s do not match. The minimum product distance approach is not suitable 
for coded schemes. For example, the minimum product distance is zero for ^ = degrees, while 
this rotation angle is optimal in terms of outage probability. 

The outage probabilities of the considered multidimensional constellations are shown in Fig. 
|8(b)[ This confirms that constellation expansion together with the optimization of the precoding 
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parameter is sufficient to approach the outage probability with an i.i.d. Gaussian input alphabet 
very closely. It also shows that the constellation = S-TZ'^ with 6' e [0 ... 9] degrees represents 
the best trade-off between performance and complexity. 

2) Extension to complex constellations: All the proofs in this paper are valid for complex 
constellations. This means that also for complex constellations, the outage region is outer bounded 
by a i?-hypersphere, determined by one variable, its radius. We restrict our attention to real- valued 
precoding matrices. Consider Eq. (|2), where z is now a complex vector and P is real-valued. 
For complex symbols, this can be rewritten as 

X = P7^{z} -t- jPI{z}, (20) 

where = —1, TZ{.} and I{.} take the real and complex part respectively. This means that the 
real and imaginary part of the complex vector are each precoded by the same matrix P. 

Assume that a transmission rate i? = 1.8 bpcu is aimed. Initially, we take the constellation 
VLz = (Rc = 0.9), which can be build as the Cartesian product of two 4-QAM constel- 

lations (16-C^=4-QAMx4-QAM). As for real-valued constellations, the rotation angle 9 can be 



optimized, see Fig. |9(a)[ The gap to the outage probability corresponding to an i.i.d. Gaussian 
input alphabet can be closed by a constellation expansion and a new optimization of the rotation 
angle. The same strategy as for real- valued constellations could be applied by only adding one bit 
in the multidimensional constellation, which would extend fi^ = to ^Iz = However, 

= cannot be written as the Cartesian product of two constellations and is therefore 

less convenient to generate (32 points would have to be placed properly in a 4-dimensional 
space). For simplicity, the constellation expansion is done by adding one bit per component 
which extends = IQ-C^ (=4-QAMx 4-QAM, = 0.9) to fi^ = 6A-C^ (8-QAMx8-QAM |], 
Rc = 0.6). The optimization of 9 and the optimized outage probabilities are shown in Fig. |9] 

Note that, for = the profile of the rotation angle 9, and so the optimum rotation 

angle, is the same as for = 4-7^^ and R = 0.9 bpcu. This can be explained as follows. 
When Tl{z} and I{z} are drawn from the same real-valued constellation then z belongs to 
a constellation = ^2 + j^z- Consequently, x belongs to a constellation = + j^x, 
where is obtained by applying the precoding matrix P to the constellation "^z- From the 

*The 8-QAM constellations have the same forni as in Fig. |6(a)| 
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Fig. 9. The optimization of 9 and tlie optimized outage probabilities are shiown whien complex symbols are transmitted. The 
transmitted rate is R — 1.8 bpcu. 



chain rule of mutual information flS], we obtain 

/n. (7) = (I) . (21) 

Hence, the precoding matrix P that is optimum for a real- valued constellation \E'^ and rate R is 
also optimum for a complex constellation fi^ = ^ + j'^z and rate 2R. The corresponding SNR 
for the complex constellation is 3 dB higher than for the real- valued constellation. For = 16- 
and R = 1.8 bpcu, the corresponding real- valued constellation is = 4-7^^. Therefore, the 
profiles in Figs. |9(a)] and |7] are the same for both constellations, except for an upward translation 
of 3 dB. 



We also tested the performance for higher spectral efficiencies. Consider for example a system 
requiring a spectral efficiency of R = 3.6 bpcu. Here, the same techniques can be used. First, 
= 256-C2 (16-QAMxl6-QAM, R^ = 0.9) is optimized, followed by = 1024-C2 (32- 
QAMx32-QAM0, Rc = 0.72). The results are given in Fig. [TOl The same observations hold as 
for the previous numerical results. 

B. Numerical results for B = 3 

When B = 3, O = 01 (see Eq. ([8])), and the optimization criterion for the upper bound on the 
outage probability is to find 9i so that I^^(9i, BR) is minimized. Next, a constellation expansion 

^The well known cross 32-QAM constellations are used. 
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is performed to further minimize the upper bound as well as the lower bound on the outage 
probability. 

We aim to transmit R = 0.9 bpcu. We construct constellations fi^ that are invariant to a cyclic 
shift of the components of the constellation points. Hence, the constellations are invariant to a 
rotation of 27r/3 with respect to the bisector (1,1,1) (this can be verified by evaluating P in 
Eq. ([8]) for this rotation angle). An example of such constellations are the Cartesian products of 
three identical one-dimensional constellations, such as 8-7^'^ =(BPSK)'^ (constellation points are 
comers of a cube) and 64-7?.'^ = (4-PAM)^ (constellation points are the comers of four nested 
cubes). More sophisticated constellations can also be considered. For example, 16-7^'^ can be 
constructed by considering 5 sets of three constellation points {(%, bj, Cj), (bj, Cj, aj), (c-,, aj, bj)} 
y for j = 1, 5, and adding a constellation point located on the bisector To build the 
actual constellations, a few design parameters have to be specified, such as the distances between 
the planes perpendicular to the bisector, the radii of the circles containing three points of the 
constellation in the planes, how many points on the bisector are taken and how many points in 
groups of three are taken. These design parameters will impact on the final performance. 

Because it is not the main topic of the paper, we will not elaborate on many different designs 

*The three constellation points are in a plane perpendicular to the bisector. 
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Fig. 11. The optimization of 61 and the optimized outage probabilities are shown for 5 = 3. The transmitted rate is R — 0.9 
bpcu. 



for the multidimensional constellations. We compare the performance of 8-71^ {Rc = 0.9), 16- 
7^^ (Rc = 0.675) and 64-7^^ (Rc = 0.45) in Fig. [Ill Note that the symmetry point in Fig. 
ll(a)| is not 9i = 45 degrees, as for B = 2, but it is 60 degrees. The S-TZ^ constellation 
is equal to the Cartesian product of three BPSK constellations and the 64-7^"^ constellation is 
the Cartesian product of three 4-PAM constellations. For the 16-7^'^ constellation, we chose 
to take a circle of three points in a plane containing the origin perpendicular to the bisector, 
and two circles with 6 points, each in a plane next to the first plane, perpendicular to the 
bisector. Finally the origin is also chosen as a constellation point. The rounded coordinates of 
the points of the constellation are {(0.0, 1.0, -1.0), (1.9, 0.12, 0.12), (1.3, -0.5, 1.3), (0.5 - 
1.3 — 1.3), (—1.9, —0.1, —0.1), (0, 0, 0)} as well as the cyclic shifts of these coordinates. 

C. Practical implications 

The motivation for this work was to solve the theoretical problem of determining the optimal 
precoding matrix P in terms of the outage probability for a given combination {Rc, ^Iz} satisfying 
a target spectral efficiency R. Besides the theoretical relevance, the numerical results suggest a 
practical relevance because the impact of this optimization on the outage probability is significant 
in the case of large coding rates Rc 1^. This case might be preferred upon the case of small 
coding rates (and thus large constellation sizes for a given spectral efficiency R) because 

'of course, when the coding rate decreases, the impact of the error-correcting code on the error rate performance increases 
yielding a decreased impact of the precoder on the latter 



January 18, 2013 



DRAFT 



25 




C I , , I , , I , , I , , I , , i_ 

7 10 13 16 19 22 

Eb/No [dB] 

Fig. 12. The optimized outage probability can be approached very closely by the WER of a practical error-correcting code. 
The block length is iV = 5000, but the results are valid for all block lengths (tested for iV = 2000 and N = 10000). 

large constellation sizes complicates the system design (requiring a joint optimization of the 
constellation labeling and error-correcting code which is not trivial) and is sometimes even not 
allowed (e.g., < 64 in 4G systems). 

The consequences of our work on practical code design are as follows. The WER of a practical 
system is lower bounded by the outage probability (the lower bound is achievable). To minimize 
the WER, first its lower bound must be minimized. Therefore, the multidimensional constellation 
^Iz and the rotation angle interval for the practical code should be taken as obtained in this work. 
Next, the labelling, the rotation angle within the rotation angle interval obtained in this paper, 
and the error-correcting code must be determined. The last three optimizations are the topics of 
another work [6|, but it is important to understand that these optimizations are based on what 
is presented in this paper. In Fig. [I2l we show for B = 2 and R = 0.9 bpcu that the optimized 
outage probability can be approached very closely by the WER of a practical system. 

VIIL Conclusions 

We have studied the effect of linear precoding on the outage probability of block fading 
channels. We have analyzed the outage boundaries in the fading space and established outer and 
inner boundaries with simple shapes which yield an easy optimization of the outage probability 
for a discrete constellation, for an arbitrary number of blocks in the fading channel, real- valued 
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or complex constellations, low or high spectral efficiency. The combination of a constellation 
expansion and an optimized precoding matrix, has shown to be sufficient to closely approach 
the outage probability corresponding to an i.i.d. Gaussian input alphabet. With this work, the 
practical code performance can be optimized by admitting the parameters obtained here. 
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A. Symmetry conditions for multidimensional constellations 

The points {0:6,0, b = 1, . . . , -B} correspond to the case that all fading gains are zero, except 
one, whose value is the scaling factor of the projection of the multidimensional constellation on 
the b-th coordinate axis Xb, so that the mutual information between X and Y is equal to the 
spectral efficiency BR. In other words, if the projection of the multidimensional constellation on 
each coordinate axis yields the same set of points, then the magnitudes of the points {0.1,^0, b = 
1, . . . ,B} are equal. 

First, we restrict our attention to the case that B = 2. Consider the constellation point z^'^ = 
{ui\u2^) G Qz- The projection of the multidimensional constellation on each coordinate axis 
yields the same set of points if for each point z^^\ the points z*^-^^ = {u^^\u2^) and z*^^^ = 
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(«i ) exist, i,j,qe[l,... 



2""]; j,q 7^ so that 



cos{e)u? - sm{9)u^^ 




In other words, = X2'' and = x^i \ where x^'^x'^-'^ and x^''^ are the corresponding points 
of z«,z(j) and z^'?) in n^. It can be easily verified that this is always fulfilled if 
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or in other words, the constellation is invariant under a rotation of 7r/2, which is obtained after a 
cyclic shift and a reflection, which proves what was claimed. An example of such a constellation 



is the constellation shown in Fig. 2(a) 



Now consider the case that B > 2. Consider the i?-dimensional constellation that contains 
M points. When z belongs to Q^, then also z*^^), . . . , tS^^^^ belong to f]^, where tS^^ is obtained 
from z by a 6-fold upward cyclic shift of the components of z: z^^^ = C^z, where C is obtained 
as a cyclic shift to the right of the columns of the B x B identity matrix. Note that the number 
of constellation points does not need to be a multiple of i?: a subset of the constellation may 
consist of an arbitrary number of constellation points of the type [z, z, . . . , zY which remain 
invariant under a cyclic shift. 

Consider an orthogonal circulant B x B precoding matrix P. Therefore, P = CPC^ (a 
circulant matrix remains the same when applying a left cyclic shift to the columns and an 
upward cyclic shift to the rows). The transformation of tS^^ is 

Pz(^) = PC\ = {CPC^)C\ = CPC^-^2 = ... = C''Pz = C^x = x^''), (22) 

where we exploit that C is an orthogonal matrix. 

Consider the matrix (x, x.^^\ . . . , x*^^"^^) . As the {i + l)-th row is obtained as a cyclic shift to 
the left of the i-th row, the set of components in a row is the same for each row. A constellation 
point in ^1^ of the type {z, z, . . . , z)^ is transformed into a constellation point in ^l^ of the type 
(x, X, . . . , x)^. We conclude that the projection of the constellation il^ on any of the coordinate 
axes yields the same set of points. 

B. Proof of Lemma [7] 

The mutual information J(q;,7,P = J) is equal to ^^/(Y; X|q:, 7). This can be split in the 
sum of B terms through the chain rule [13 : 

J (a, 7, P = J) = 4 {I {Yi; X|ai, 7) + J {Y2; Xl^i, ai, «2, 7) + • • • + ^ (Yb; ^\Yi, . . . , l^-i, a, 7)) 

JD 
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The mutual information 

/ (n; xiFi, . . . , n-1, cx,^) = H (nin, . . . , n-i, 7) - ^ (n|x, . . . , n-i, 7) 

= H (nii^i, . . . , n-i, a, 7) - ^ {y^x,, x^,, y,,..., y,.,, a, 7) 

= H {Y,\Y,, . . . , n^i, a, 7) - i/ {Y,\X,, Fi, . . . , a, 7) 
= /(n;X,|ri,...,n-i,a,7), 
where X^;, = [Xi, . . . , X^+i, . . . , X^]. Because 

H iY,\Y,, ...,Y,_,,(x,^)<H {Y,\a, 7) (23) 

and 

H (n|X,, Y,,..., n-1, a, 7) = (niXfc, a, 7) , (24) 

it is clear that 

/ (X,; nin, . . . , «, 7) < / (X,; 7) , (25) 

where the upper bound equals Ig^ («^7, P). 

C. Proof of Props. |2] and \3\ 

Consider a function F{x) that is concave for a; > 0. Hence, for arbitrary L, 

\/=i / 1=1 

for J2i=iPi = 1, Pi ^ ^rid 7i > for / = 1, . . . , L, by Jensen's inequality. In addition, assume 
that F(0) = 0. 

We construct the function J2b=i ^ilb) where Ylb=i lb = C, ■yb > and 7^ = ■yaf for 
6 = 1, . . . , i?. Hence, a is on the surface of a hypersphere with squared radius C/7. 
From Eq. (|26l ) with L = B, we obtain 

E^(7.)=5(^;^i:i^(7.)) <fii^ (^1^^') ^^'^ (^) 

The maximum value B.F (^) is achieved for 7^ = ^, b = 1, . . . , B. 
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Further, using Eq. (|26] ) with L = 2, 



Summing over b yields 



I > (C) + (0) = (C) 



b=l 6=1 

The minimum value F{C) is achieved when, for given / G{l,...,i?}, 7z = C and 76 = 0, for 
67^ /. 



D. Proo/ 0/ Prop. |?] 



rrr = 7|Q;p. In this proof, we assume real constellations 



The instantaneous SNR is = =771— j 
for simplicity. The proof can be straightforwardly extended to complex extensions. We denote 
the real part of the complex noise w by w'. 

Consider the mutual information / (ck, 7, P) of the constellation Qt, given the fading gains a 
(Eq. (flOl)). This expression can be reformulated in terms of 7: 



/K7,P) = ^-^$^Ey|. 



B B 



log2 XI (^^(y' « ■ - c^^(y, a ■ x'))] 



where Ey|x can be replaced by an expectation over the noise, E^', w[ ~ A/'(0, 1/(27)). The 
argument of the exponential functions can be simplified, so that 

B 



m 2 

/ (a, 7, P) = — - -g- ^ 



10& E 



exp 



-7rf'(a -x, a ■x') + 5^(27^^(6)) 



6=1 



where /(6) = ab{xb — x^). This expression can be further simplified by approximating the 
exponential functions and logarithms by their respective Taylor series for small 70^, V6 = 
1, . . . , B. Next, the expectation of the expression over the random vector w' can be replaced 
by an expectation over 7w', where 7w' ~ A/'(0, |I) (I is the identity matrix). Therefore, we 
can drop all terms that are proportional to E^^' [jw'f] and replace E^^' [(7^^^,)^] in all terms 



proportional to E^^/ [(T^t,) ] by f • Now, after some calculus,we obtain 



B 



T( u\ 7 E6 "6Var(X;,) ^ 
/ (a, 7, P) = %^ZZ777\ + 2^ 



51og(2) 



(27) 



6=1 
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where Var(Xb) is the variance of the 6-th component of the points of constellation Q^. The 
validity of (ITTI) for small instantaneous SNR has been verified numerically. As the projection 
of on either coordinate axis yields the same set of points, this variance is independent of b. 
Hence, for small instantaneous SNR, the mutual information remains constant for the set of fading 
points where J2b=i ^1 constant and the outage boundary coincides with a i?-hypersphere. As a 
consequence, and V^tte are equal and by their definition, it is clear that for low instantaneous 
SNR, the outage region coincides with the 5-hypersphere J2b=i '^6 = '^o = \fBa^. 

E. Proof of Prop. [5| 

Recall that I (a.,'y,P) = -^/(T; Y|a, 7). For high instantaneous SNR, it was conjectured in 
[|27l . [|26l that maximizing (minimizing) the minimal Euclidean distance dmm of a constellation 
maximizes (minimizes) the mutual information on a Gaussian channel. This was based on tight 
upper and lower bounds on the mutual information, established in ETl . [|26ll (generalizing [jlD ). 
More specifically, for dmmioi) > 0, 

where dmmict) is the minimal distance of the constellation Clt, and where Ki,K2, and K^, are 
constants. When d'^^^{c3L) = 0, the constellation Qt effectively counts less than 2™ points and the 
mutual information /(T; Y|q;, 7) will converge to a value strictly smaller than m for increasing 
SNR. We now prove that, for cki, 0.2, satisfying d'^:^^{a2) > ^^^^(ai), there exists an SNR-value 
7 above which it is impossible to have 

J(T;Y|ai,7) >/(T;Y|a2,7). (29) 

This is obvious for dj^j„(ai) = 0. Let us check whether (|29|) is invalid when d'^^^(a.i) > by 
using contradiction. More specifically, let us assume that (|29|) is true. Combining with (|28T ). it 
follows that 

e-'^min("i)7/4 / \ e~'^^'"("2)'>'/^ 

^ C?min(ai)V7 V ^ <„("l)7/ ^ ^C?mi„("2)v^' 
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yielding 

I (C?min("2) - C?min("l)) 

The right-hand side of (|301 ) is a constant close to ^^^"'p'^"^^ for large SNR. Clearly, there exists 



a large SNR value for which (pOt . and thus (1291 ) is not true, confirming the conjecture made in 

flan, m. 

Outer bound: 

Consider two points x*^*^ = {u(\ . . . , u"^) and x*^-^^ = {uf \ . . . , m^^) from the constellation Vt^. 
The corresponding points t^*^ and t*^-'^ from the faded constellation Vtt have a squared Euclidean 
distance given by 

B 



Y.ai(v!i^-v!i^) . (31) 



2 



Let us denote by h* the value of 6 G {1, . . . , 5} for which ^m^*-* — m^"*^ j is minimum. For the 
B-hypersphere Xlf=i "^6 = '^o' ^e obtain 

|t«-t(^Y>«o(4^-«i'0'' ^^^^ 

where equality in Eq. (|32] ) is achieved when a;,* = cto- When 6* is not unique, equality in Eq. 
(|32|) holds when a?,* = ao holds for any h* that minimizes — j . Hence, the minimum 
distance dmm{oi.) for the constellation Vtt is given by 



dmmioi) = cto min min 
i,ie{i,...,M} be{i,...,B} 



(33) 



As the constellation is such that its projection on either coordinate axis yields the same set 
of points, the minimum distance (/^^^(q;) on the i?-hypersphere is achieved for either fading gain 
equal to ao (and the remaining B — 1 fading gains equal to 0). Hence, according to relation 
between minimal distance and mutual information, I(ol,^,P) with a on a i?-hypersphere is 
minimized on the coordinate axes at high SNR. Taking so that I(a.,'j,P) = R on the 
coordinate axes yields / (ck, ^, P) > R on the other points of the i?-hypersphere. 
Inner bound: 
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We must prove that the minimal distance of Q,t is maximized in a.^, or 



argmax mm 



iiin > 



au u 



For OL = cte, we obtain J2b=i 



at u 



u 



6=1 
2 



u 



0^1 EL 



u 



(i) 



U 



(34) 



which is mini- 



mized to a^cP when x*^*^ and x*^^^ are points at minimum distance d in the constellation 



< 



For any a with |ap = Ba^, a ^ a^, we will prove that minj j ^^^^ 
Ofgd^, which is what is claimed. We first restrict our attention to the case B > 2. We con- 
sider B pairs | ^x*^"\ x'^"-* j , n = 0, . . . , i? — l| of points in and their corresponding points 
|(^t('^),t'^''^) ,n = 0,...,5-l| in Qt, such that x^'^) = T"x and x'^"^ = T"x', with T 
denoting the matrix operator that causes a downward cyclic shift. We select x and x' such that 
|x — x'l = d; obviously, jx^"^ — x'*-"^! = d for n = 0, . . . , B — 1 . Because J2b=i (^^''^^^ 

^/(n)|2 



,("■) 



'(«) 



)z 
, we obtain that, for any ct with jcKp = Ba"^, the average distance E„|t(")- 




2 / W '(n) 



1 ^ 

6=1 



B-1 
n=0 



(n) 



'(n) 



d2 



j2 2 

d az. 



(35) 



7 



'(n) 



Hence, min„ Y.b=i "fe ( 
we consider 4 pairs 

and x'*^"^ = T"x', with T denoting the matrix operator representing a counter clockwise rotation 



< d'^a^. For B = 2, the proof has to be slightly modified: 
) , n = 0, . . . , 5 - l| of points in fi^., such that x(") = T"x 



of 7r/2. We select x and x' such that |x — x'| = d; obviously, |x^") 
For any a. with |q;|^ = 2al, we have, considering that = — /, 

/ 



I 



d for n = 0, 



3. 



n=0 \b=0 



au U 



,(«) 



u 



'(n) 



IE 



6=0 



E 

n=0 



(n) 



U 



'(n) 



2d2 



j2 2 



(36) 



The remainder of the proof follows the same lines as for B > 2. 



January 18, 2013 



DRAFT 



33 



F. Proof of full diversity 

Let us find an expression for a^- Tlierefore, we consider that all fading gains are zero, but 
one, so that the rate achieved by Vtt (which is now a scaling of iSp) over that single slot must 
be at least RB, or Is^{al'^, P) = BR, so that 

2 ^s:iBR,P) 
«o = — , (37) 

7 

Now, as the outage region Vo is outer bounded by the hypersphere with radius cto, we have 

Pout(7, P, R) < Pr(a? + al + ... + al<al). (38) 

Because the fading gains {a^, b = are i.i.d. Rayleigh, the cumulative distribution 

function of the chi-square distribution with 2B degrees of freedom is fl28l 

B-l 

Pr(«? + al + . . . + al < x) = 1 - e-'^y^lJ- 

fc=0 

The right hand side can be manipulated as follows: 

B-l u X x*= V^oo rf.B Y^oo 

1 ^ ^ g - Z^fc=0 fc! ^ l^k=Bld ^ Z^fc=0 {B+k)\ B 

fe=o Z^m=0 m! Z^m=0 m! 

so that from Eqs. dSV]) and (|38]) 

Pout(7,^,^)<ao''«^- (39) 

Hence, with increasing 7 the outage probability goes to zero not slower than 7"^, so the diversity 
order is at least B. As there are only B i.i.d. fading gains involved in the transmission of a 
codeword, the diversity order is exactly B. 
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